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We describe a recent progress in finding solutions to three-particle evolution equa- 
tions at leading order in the QCD coupling constant for multiparton correlation 
functions based on the integrability of corresponding interaction Hamiltonians. 

1 Multiparton correlations in hard scattering. 

A typical lepton-hadron cross section at high momentum transfer Q is given 
by a series in the latter 



where r stands for the twist of contributing operators which parametrize the 
physics at soft scales. The first term in the expansion corresponds to the 
Feynman model picture. Consequent terms, standing for multi-parton corre- 
lations in hadron, reveal QCD dynamics not encoded into the conventional 
parton densities. They manifest genuine quantum mechanical effects being an 
interference of hadron wave functions with different number of particles, see 
Fig. 0. Twist-three correlations are unique among other higher-twist effects 
by their feature to appear as a leading effect in certain spin asymmetries. The 
most familiar examples are the transverse spin structure function 52 measured 
in deep inelastic scattering. If a cross section is measured in a rather wide 
range of hard momentum Q, the logarithmic scaling violation in the functions 
/ in Eq. (0) becomes an issue since hadron substructure is seen by a probe 
with different resolutions. In QCD it is formalized by means of renormaliza- 
tion group equations with kernels ■wiucb are derived by standard methods of 
QCD perturbation theory (see e.g. BB'El). The specifics of twist-three sector 
as compared to the twist-two case is that the afore mentioned structure func- 
tions being defined as Fourier transform (FT) of two-particle hadron-hadron 
matrix elements / ~ FT {h\(j)i(j)2\h) receive contributions from three-parton 
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Figure 1: Hadronic structure as a series in a number of Fock components participating in 
a hard rescattering. The first term on the r.h.s. of the equality stands for leading twist-two 
contributions. The second one (plus its complex conjugate), to the twist-three effects, etc. 



quark-gluon correlators ui3'u plus a kinematical piece from twist-two operators 
going under the name of Wandzura-Wilczek part. The original operators, not 
possessing a definite twist, mix with two- as well as three-particle operators 
and thus evolution equation will have an extremely complicated form. Resolv- 
ing this complication in favour of an independent study of the renormalization 
group evolutions of separate twist components one finds an autonomous equa- 
tions for two- and three-particle sector, respectively. The twist-two sector is 
simple and is well-studied in the literature. It is the latter which is the subject 
of this presentation. In spite of simplifications due to reduction alluded to 
before, anyway, one ends up with a very complicated problem of working out 
the mixing of three-particle local operators. 



2 Leading order evolution. 

As we have stated above reducing the twist-three two-particle operators to the 
three-particle ones 0i02</'3 allows to easily find corresponding evolution equa- 

ton since the latter operators fall into a class of the so-called quasi-partonic 
ones and thus have a bunch of remarkable properties which simplifies their 
renormalization properties. Namely, at leading order in the QCD coupling 
constant the total kernel if 123 reduces to the sum of conventional pair-wise 
twist-two non-forward kernels Kab which includes the momentum conservation 
delta function 5{ka + kb — k'^ — k'^^. Thus, the leading order evolution equation 
is of the form 

dlnQ^ ^^^^ = ~^-^i23 * -F123, with K123 ^ K12 + K23 + K13. (2) 

Here F123 = F{xi, X2, X3) and * = J dxidx2dx3S {xi + X2 + X3 — rj) and vari- 
able 77 which encodes the skewedness of the matrix element: rj = Q for the 
forward scattering and 77 = 1 for totally exclusive kinematics. The solution to 
Eq. (S) is expected to be of the form 



Fi23(Q^)-E'^wi23 PtSt m.}{Qi))), (3) 

{a} ^"^^^ 
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where S{a} the eigenvalues of evolution kernels (and /(c) is an extracted 
colour factor) and tf'{Q,}(x, x') are the corresponding eigenfunctions paramet- 
rized by a set quantum numbers {a}. As usual /3o = — ^-Ca is the 
leading term of the QCD /3-function and {{^{a}{Q'o))) stand for reduced matrix 
elements of local operators at a low normalization point Qq. Therefore the 
evolution equation (||) can be reduced to a stationary Schrodinger equation 

-f^l23 * ^{a} 123 = f(c)^{a}'I^{a} 123- (4) 

Simplifications in diagonalization of this equation occur due to use of QCD 
symmetries which we now use in turn. 

3 Conformal symmetry. 

The classical QCD Lagrangian enjoys the property of invariance under the 15- 
parameter group of conformal transformations which consists of the Lorentz 
group, and A4piy, and generators of dilatation T> and special conformal K,^ 
transformations El. At quantum level the trace of the energy momentum ten- 
sor develops an anomaly 0^^ ^ and the generators V and /C^ cease to be 
symmetries of the theory. The invariance of Lagrangian under a given transfor- 
mation imposes severe restrictions on the form of eigenfunctions of composite 
operators entering the evolution. For illustration purposes let us consider the 
situation with two-particle operators since they enter as a building block in 
(a6)-subchannel for our general consideration of three-particle problem. It was 
known for a long time that these are conformal operators □ which diagonalize 
the leading order renormalization group equation for local operators with total 
derivatives, e.g. (jibida + db)X{da — dbY+cpa, where the -I— subscript stands for 
projection with a light-like vector n^, i.e. u+ = v^n^. Let us demonstrate it 
explicitly. The conformal operators have the form 

0,k - Mida + idbf+P'-"""^^ {{da - db) + /{da + db)+) (j^a, (5) 

where pj"^*"""^ 

a,rG Jacobi polynomicils with — c^^-I-Sq — 1, and da and Sa being 
dimension and spin of a field (pa ■ Since we deal with operators on the light-cone 
it is enough to consider a coUinear sub-group of the conformal group, which 

consists of projections 7^+ = V^n^, M |- = M^^n*^n^., V and K- = K.^n*^, 

forming an so(2, 1) algebra, [J'^,J^]- = ±J7=^, [J+,J^]- = -2J^, defined 
by generators of the momentum — i'P+j special conformal transformation 
J'^ ~ -^IC- and a particular combination of dilatation and angular momentum 
J3 ^ ^(^x> + M-+). The Casimir operator is = J^{J^ - 1) - J^J~- 
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Simple calculation shows that the operators (||) transform covariantly under 
dilatation and special conformal transformation, namely, 

t[Oji,V]^^-{l + da + db)Oji, t[Oji,IC-]-=ia{j,l)Oji-i, (6) 

with cLjkil) — Sjk ■ a{j, I) and a(j, I) — 2{j — + I + Va + i^b)- From these 
equations one easily sees the meaning for J'^ to be step-up and -down operators 
as V+ adds one unit of spin to the operator in a conformal tower while /C_ 
reduces it by one. These operators possess conformal spin Jab ^ j + ^{>^a + 
Vb + 2), [j'ib^ Ojk]- = Jab{Jab - '^)Ojk, spin I + Sa + Sb and scale dimension 
I + da + db- Thus conformal operators form an irreducible infinite dimensional 
representation of the collinear conformal group spanned by bilocal operators. 
The renormalization group equation for conformal operators have the form 

j^0,--±,,.0,,, (7) 

with triangular matrix of anomalous dimensions ^jk, j > k. To find a con- 
straint on the form of the anomalous dimension matrix of Ojk stemmed from 
the special conformal symmetry, we use Ward identities for V and /C_ trans- 
formations. Since the dilatation WI is equivalent to the renormalization group 
equation we get an equation on jj^ . Explicit analysis gives Q 

Ijj' a-fkil) - ajf (I) -fj'k = 0, (8) 

at leading order in coupling constant, 7 ^ 0{as). Since the matrix a in ^) is 
diagonal it forces 7 to be diagonal as well at lowest order in strong constant, 
Ijk — Ij^jk- Beyond this order conformal symmetry is violated and the diag- 
onal matrix a is promoted to a non-diagonal oneB ajk — > a-jk +7*^ + 2^ 6, with 
new objects which are special conformal anomaly 7"^ appearing from the renor- 
malization of the trace anomaly and a conformal operator Oji ^ dx x^O^^^ oc 
'Y^'k^ol'jk^ki, and an a^-independent matrix h. These extra terms induce the 
non-diagonal part of the anomalous dimensions beyond leading order. 

Therefore, the net product of analysis is that due to preservation of confor- 
mal invariance for the leading order anomalous dimensions (pair-wise kernels 
of the previous section H) one immediately concludes that the former are the 
function of eigenvalues of the Casimir operator JT"^ . Thus the interaction ker- 
nels can be expressed as functions of Casimir operator in a given subchannel, 
e.g. Kab = KJl,). 
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Let us give an explicit example of a basis, convenient for present studies. 
Choose a space V = {9^\k = 0, 1, . . . , 00} spanned by elements 9'' = d^(j)/T{k+ 
ly+l). In this representation [J^'^, x(6')]- = J^'^x{9) the generators take the 
following form 

J+ ^{„ + l)9 + 9^dg, J-=dg, ^ ^{ly + 1) + 9dg, (9) 

where dg = and the quadratic Casimir operator reads = J^{J^ ~ 1) — 
J~^J~. The advantage of this basis lies in the fact that conformal operator has 
the form of translation, so that the highest weight vector condition is easier to 
solve. For a multi- variable function x{9) = x(^i, 6*2, • • ■ , 9n) the operators are 
defined by the sum of single particle ones as J^'^ = ^e^'^ ^^'^ they obey 

the usual commutation rules [J^, J"*"]- = ^J"^, [J^ , J^]- — ~2J^. Obviously, 
a single particle state 9i is an eigenstate of the Casimir operator J'^9a — 
ja9a, with = — 1). The eigenstates of two-particle Casimir operator 
in (a5)-subchannel are {Jab)''9ab (here and throughout 9ab = 9a — 9^) 
and coincide with the bilinear conformal operator (^) and possess the same 
eigenvalues. Now, the three-particle basis efficient for present applications has 
to be chosen so that it diagonalizes the total three-particle Casimir operator 
= J\2 + J\i + Jiz — J2e=i3j ^-nd one in a subchannel J^^, say a — 1 and 
b = 2, 

J^'Pj-.j = {J +\{vi+i^2 + + 3)) (J + + 1^2 + 1^3 + 1)) Vj,j, 
JliVj.j = {] + 1^2)) (.? + \{vi+iy2 + 2)) Vj.j. (10) 

The solution is given in terms of hypergeometric function T^jy (6*1, 6*216*3) oc 
9i29'-^2FiU - J, J + i/i + 1, 2j + y^ + u2 + 2\9) with 9 = 012/^32. 

Since the physical anomalous dimensions have to be real we have to define 
an appropriate scalar product resulting in selfadjoint Hamiltonian, namely 

{X{9)\X{9))^ /, \{'^^{l~9Ar-\{9)x{9). (11) 

and 9 — 9*. Then it can be seen that the Casimir operator and, as a con- 
sequence, the Hamiltonian are selfadjoint operators w.r.t. such defined inner 
product (J^)^ = J^, /i^ = /i, and therefore Imf = 0. 

Obviously, to solve the three particle problem (^) one has to find an extra 
integral of motion in addition to the one provided by the conformal symmetry. 
As we show in the next sections it turns out that this is the case of antiquark- 
gluon-quark and three-gluon systems in some limits. 
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4 Antiquark-gluon-quark correlators. 

The antiquark-giuon-quark correlations contribute to a number of single spin 
asyrnnietiesuparujLenter as a three-particle part in the twist-3 (skewed) parton 
distributions BI13'0. For instance, the transverse spin structure function g2 
admits the following three-parton piece 

gf-^[x)^ J dxidx3J{xi,X3)Y{xi,X3), (12) 

where JT^ is a differential operator acting on the momentum fractions Xi of 
a two-argument function Y(xi,X3) which is defined as a light-cone Fourier 
transform of a hadron matrix element of a C-even combination (p|iS+(ki, K3) + 
iS~(— K3, — Ki)|p) of nonlocal operators S^{ki,K3) = igip{K3n)j+[iG^^(0) ± 
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75G^_i_(0)]'(/'(Kin). This will be the case of study in this sectionEj^ 
particular we elaborate on operator 5+, since is related tn it by charge 
conjugation. Similar considerations apply to chiral odd sector E3£3'Eil. 

In view of our discussion in Section |^ the total Kqgq evolution kernel 
splits into a sum Kqg + Kgq + Kqq ~ j(3o, where the /?- function term is due 
to presence of the coupling constant in the definition of the operators S. The 
leading twist QCD interaction kernels Kab with npn-zero momentum transfer 
in i-channel are known for a long time, see e.g. □. One of the complications 
to solve the three-particle system is due to non-trivial colour structures of the 
latter, e.g. Kqq depends on Cf — '^Ca, while K(^q_q)g on both Cf — '^Ca and Cf- 
Obviously, the kernel simplifies drastically in multicolour limit N —>■ 00^ since 
the interaction of end point quarks effectively vanishes, Cf — \Ca ^ 0{N^'^). 
In this case the total interaction kernel takes a form Kqgq = ^NcTi where 

^ ha il) + t^gg ih) - h with hab{S) ^ ^{Jab + S) + ^P{Jab ^ S) - 2^{1), 

(13) 

and = J {J — 1). Thus we have reduced the original problem to an open 
chain of three particles on a line with conformal invariant interaction of neigh- 
bors (quarks and gluons) and no interaction of end points (quarks). This 
observation turns out to be fruitful since one can find the pair-wise Hamilto- 
nians ( |l3|) among the series of integrable ones for inhornDgeneous spin chains. 
It can be easily generated from the Yang-Baxter bundlellj, which is a solution 
to the Yang-Baxter equation EZl, 

Rab{\) = f{v, A^Pab-^^ ^ f 14 

T{Jab-\)T{v + l + \) 

as follows hab{6) = Rab{—S)R'^^^{—5), provided the function / is chosen like this 
/(i/, A) — r(i^A) r^tz+i-A) • total H. can be produced from the open spin 
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chain transfer matrixESl ih(A) = trf, Tb{\)Tj^ (^'^)L -2 constructed out of i?'s, 
Tft(A) = i?a^b(A — Si)Ra2bWRa3b{^ ^ "^s) with the spins of spaces chosen 
Ukc i^i = 1/3 — 1, 1^2 — using the formula Ti. = ^^lntf,(A). Identifying 
imhomogeneities as 5i = | and S3 — ^ wc immediately find Eq. (|l3|) . This tells 
us that the antiquark-gluon-quark system is exactly integrable in multicolour 
limit. This fact implies that there exist a family of commuting integrals of 
motion. The first of them is provided by the conformal symmetry of interaction 
and is given by Casimir operator j^. Thus, we have to find only one 'hidden' 
conserved charge. This follows from the standard formalism of integrable spin 
chain models and forces us to calculate transfer matrix with auxiliary space 
being two-dimensional. WiliiiihDmogeneity parameters taken as above it gives 



an expansion in rapidity AcJiijO 

t^{X) - r2(A) - (4A2 - l){X^-3l)J^ - i(4A2 - l)Q{S^,Ss) , (15) 
with constant function i7 and the 'hidden' charge 

Q {61,63) = [Jj^Jlsh - ^^iJh - 2SlJi2- (16) 

Below we will find eigenfunctions of this much more simple (compared to orig- 
inal Hamiltonian) operator Q and use them to find the eigenvalues of the 
Hamiltonian (^3|). Let us note in passing that integrable open spin chains have 
been encounteced as well in the solution of BFKL-type equation for quark- 
gluon reggeonslia. 



To find the eigenfunctions of the Hamiltonian (13) we solve the equation 
for the charge with eigenfunctions 'F expanded in the conformal basis discussed 
in Section || 

J 

QilD'^'^is^^ ^ = Y.ejr,Vj,,{e,,e2\0^), (i7) 

with Qj^ = [{J -j + 1)(J + j + 5)(j + l)3(j + + 2)3] being normal- 

ization coefficients. The main advantage of the basis is a three-diagonal form 
of the matrix elements of non-proper two-particle Casimir operators »/(i,2)3- 
Then one establishes that Eq. (|l^) is equivalent to the three-term recursion 
relation 

(2j+3)r,+i + (2j+5)r,_i+e2 (2j-f3)(2j+5) ([Q (|, i)]^.^. - qs) T, = 0, (18) 

where [Q]^ j are the diagonal elements of the charge {Vj-j\Q\Vj-j) in the basis 
Vj-j{di, 6*2 1 6*3). Polynomiality requires the expansion coefficients, Tj, to satisfy 
the boundary conditions T_i — Tj+i — 0. 
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The knowledge of eigenfunctions allows to find the energy of the antiquark- 
gluon-quark system via the formula 



(19) 

with e{j) — + 1) + ^{j + 4) — 2-0(1). Unfortunately, the exact solution to 
(H), and thus (|l|), can hardly be found. Therefore, we restrict ourselves to 
WKB approximation for large conformal spins J. pEjpr the levels behaving as 
£ cx 21n J we can find the exact^ lowest trajectory E3 

£{J) = V(J + 3) + V(J + 4) - 20(1) - i. (20) 

The remainder of the spectrum is described by the formula 000 

£{ J, qs) = 2 In J - 40(1) + 2 Re0 (| + t^s) - §, (21) 

with rjs — \\J 2qs I J"^ — 3. The conserved charge is quantized, with WKB 
quantization condition arisen from the matching of the WKB solution with 
exact ones at 'turning' points, and gives quantized values of energy via (pl|). 
It compares well with an explicit numerical diagonalization of the anomalous 
dimension matrix u as shown in Figs. |^. There is an alternative description 
of the spectrum by trajectories which behave as £ oc 4 In J at large J. The 
dependence on the conserved charge reads asymptotically f ( J, qs) — In \<ls ~ 
40(1) — I + (!?( J^^). JThfi first few expansion coefficients of energy in the series 
in Jzi was found inEjO. Corrections in have been discussed recently in 
Ref.B. 



5 Three-gluon correlators. 

Three-gluon correlators enter the g2 via the quark loop coupling shown in Fig. 
^ and might be responsible for its small- a; behaviour, 

52(2;)= J dxidx3Cggg{xi,X3)G{xi,X3). (22) 

Here C,^, = C^f^ + cf^, with C^f^ calculated in Ref.i but contact-type con- 
tributions (6) not accounted for. The function G{xi,X2) is a Fourier transform 
of non-local operators, e.g. 5^(ki, K3) = g/afcc G+i('«3")G^(0)G'^|i(Kin). 

"^For the chiral odd distribu|fc«)n similar result has been found in Refs.Ej'Ej and for twist-three 
fragmentation functions incil. 
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The analysis, equivalent to the one done for the antiquark-gluon-quark 
correlator, gives the evolution kernel Kggg which can be expressed up to cor- 
rections playing negligible role in the generation of the spectrum of anomalous 
dimensions as a sum Kggg = ^CaH + |/?o where Ti. is split into two pieces, 

H = Ho + V, where = ^-12 + /123 + V = 1^12 +1^23, (23) 

and pair-wise interactions given by 

hab -^2^, (Jab) ~2^P{1), Vab = -^Jab- (24) 

Obviously, the separation into TLq and V is not accidental, but reflect the fact 
that Tioj which describes the gluon interaction with aligned hclicities, coincides 
with the Hamiltonian of XXX closed magnet of noncompact spin — |. This 



can be easily found by means of the same R matrix (14) from the logarithmic 
derivative, Ho = ^liii6('^)i of the transfer matrix of an open spin chain 
tb{X) = trbRaibWRa2h{^)Ra3bW- As a result the system described by Hq 
contains a family of mutually commuting conserved charges which arise from 
the expansion of the auxiliary transfer matrix 

ti{X) ^2X^ + X(f -if) +iQ, with Q = §[Ji2,'^23]-- (25) 



Note that the same Hamiltonian, however nf different conformal spin, has been 
encountered in interaction of QCD reggeonsEiiEa and Brodsky-Lepage evolution 
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(«) CO 

Figure 3: Coefficient function of thrcc-gluon state in the structure function 92 • 

equation for baryonsEl in the case of the same hehcities of participating quarks. 
The coincidence of quark and gluon kernels for the same heUcity states is a meie 
consequence of A/" = 1 supersymmetry which relates corresponding kernelsE3, 
/i^,(2n + 1) ^ /ij<,(2n + 1) = h^g{2n + 1) where /i^g(J) = 21/^(7 + 2) + 2V'(J + 
3)-4V'(l) + 2(-l)'V(J + 2)-3. 

So our strategy will be, first, to diagonalize the exactly solvable interaction 
and then to consider the effect of V interaction on the spectrum. Thus in 
parallel to previous section, we diagonalize first the charge Q and find its 
eigenfunctions T, (0), now gj^ = [{j + l)i{J + j + 6)2{J - j + 1)2/ i2j + 5)]^/^ , 
from the recursion relation 

2(2j + 5) qTj - a,r,-+i + a,+iT,_i , (26) 

with coefficients ttj = [j + l)(j + 3)(J + j + 6)(J — j + 2). T's allow to find 
the energy via the formula 

go(J,g) = 2Rc ^^-f '^-^ ^^ + 3, (27) 

with the two-particle interaction e(j) ~ 2ip{j+3)—2tJj{l). The lowest trajectory 
can be found exactly 

£0(^,0) = 2V'(iJ + 3) +2V'(iJ + 2) -4V'(l) + 4. (28) 

While the rest of the spectrum is described by WKB formula 

£o{J, g) = 41nry - 6V^(1) + 2ReV' (| - iq*) , (29) 

where q* = q/rf and is valid with 0{ri^^). We use the convention 7/^ = 
(j + |(i^ + 1)) ( J + |(i^ + 1) — 1) for the eigenvalues of the quadratic Casimir 
operator of the chain which to the stated accuracy reads for gluons 77 = J + 4 + 
0{ J^^). Matching of the WKB and exact solutions to (p6| ) gives a quantization 
conditions for the charge, q*h\'q = arg T (| + iq*^ + ^ {2m + i[l — (— l)"']). 
An alternative set of trajectories behaves like f (J, q) — 2 Ing — 67/;(l) + 0(J^^). 
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The perturbation V does not affect a bullc of the spectrum except for 
the lowest few levels. To analyze the situation one conskleis Hamiltonian 
Ha —'H + aV specified by a coupling constant a and study ESO the level flow 
as a function of a. Making use of these results we design the following formulae 
which describe well the spectrum of the perturbed Hamiltonian, see Fig. H, 

£{J,m) = £o(J,0)(<5„,o + 6ms) + £o{,J,q{m))e{m - 2) - A(to). (30) 

Here 5o(J, 0) is the ground state energy ( p8|) and £o{J,q{m)) is Eq. ( p9|) with 
the q{m) trajectories deduced from the quantization condition. Finally, 

A(m) = 0.54 (5„,,o + 0.08 6„,,i + {So -d{m-2))0 {Sq - 6{m - 2)) , (31) 

is the shift-function with 6q = 0.15 and 6 = 0.01. 



6 Conclusions 



To conclude, it remains a challenge fWr^.- 
p£ pintegrable structures for light-cone t3^E3ll^ 



.^understand the appearance 
and transverse momentum 



I dynamics and unravel their striking similarity. 
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